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Abstract 



We show that a radiative modification of a recently proposed model by Altarelli and Feruglio 
with softly broken A4 symmetry leads naturally to nonvanishing Ue3 with 613 ~ 2° — 4°. 
The observed mass squared differences for solar and atmospheric neutrinos are reproduced, 
whereas the predicted solar neutrino mixing angle is brought down from the tri-bimaximal 
prediction to be in better agreement with the latest global analysis including experimental 
data from KamLAND and SNO. 



Experimental measurements on neutrino oscillations are consistent with nearly maximal 
atmospheric neutrino mixing angle (^^23 — 45°), large but less than maximal solar neutrino 
mixing angle {612 — 34°), and a small "CHOOZ" angle (^^13 < 10°). Whereas the three 
neutrino masses could be quasi-degenerate or hierarchical with normal or inverted ordering, 
these values of the mixing angles are found to be remarkably close to the conjectured tri- 
bimaximal mixing ansatz of Harrison, Perkins, and Scott (HPS) 



V 



The importance of the non-Abelian discrete symmetry group A4 in understanding why 
charged leptons may have very different masses and yet a symmetry exists for the neu- 
trino mass matrix has been discussed by Ma and collaborators in recent papers |21 El El- 
In particular, it was shown in 0] how the HPS ansatz may be realized. In an interesting 
development, Altarelli and Feruglio (AF) have proposed the simplest such model with only 
two parameters in supersymmetric as well as nonsupersymmetric cases The relevant A4 
symmetric Lagrangian, after spontaneous symmetry breaking, becomes 

+ XaVl{u/ K^){lJel'e + 2z/^Z/^) 

+ I^tI^t - VeV^ - I'^.I't - l^r^^e) + h.C. (1) 

Here A is the scale of new physics (= seesaw scale) and the Higgs content with their vacuum 
expectation values are presented in Table 1 where all fields except xt have been used in Ref. 
[Sj. Then the neutrino mass matrix takes the form 

/a + 2d/3 -d/3 -d/3 \ 
M^^ = mo -d/3 2d/3 a -d/3 , (2) 
^ -d/3 a -d/3 2d/3 / 

where a = 2xau/A, d = 2xbVs/A, rriQ = v^/A. This is exactly diagonalized |1] by the HPS 
matrix, resulting from the underlying symmetry with 
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Table 1: List of fermion and scalar fields used in this model. 



where the common mass factor ttiq has been absorbed into a and d. 

Although this model prediction of 6^^ = is consistent with the CHOOZ - Palo Verde 
upper bound, sin6'5'3 < 0.16 (6'i3 < 10°) [21 Ejj the actual value of the mixing angle is of 
considerable theoretical and experimental interest as more accurate values on the parameter 
are expected to emerge from long baseline and future reactor experiments such as Double 
CHOOZ, Triple CHOOZ, Noi^a and others jH]. Thus, it is worthwhile to study some modi- 
fication of the AF model which may predict ^13 7^ 0, which is necessary for CP violation in 
neutrino oscillations. The second observation is that the tri-bimaximal mixing matrix pre- 
dicts tan^ = 0-5 corresponding to 6'°2 = 35.3° whereas a recent global analysis including 
the KamLAND and the latest SNO data gives (tan^ Oi2)expt. = 0.45 ± 0.05 corresponding to 
{9i2)expt. = 33.8° ± 1.5° f^. Within la the tri-bimaximal mixing prediction just touches the 
upper limit of the solar neutrino mixing angle and it is desirable to have a model prediction 
where the mixing angle is in accord with the central value obtained from global analysis. 

Further, while fitting the available neutrino data the following relations between \a\ and 
\d\ have been found useful. In the AF parametrisation both Aml^^ and AttIq have been 
fitted with an ansatz, 

|c/| = -2|a|(l-2/?)cos0. (5) 

where = arg((i) — arg(a) and R = Amg/Am^^^ and this relation manifests a moderate fine 
tuning of the two parameters. On the other hand, one may assume a simpler relationship, 

\d\ = —2\a\ COS0, (6) 

in which case |mi| = \m2\ at the seesaw scale. In such a scenario, while the value of Am^^^ is 
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fitted, ArriQ is expected to be generated by radiative corrections. However, such corrections 
are negligible for small neutrino masses |mj| << 0.3 eV in the Standard Model (SM) or in 
the Minimal Supersymmetric Standard Model (MSSM). New particles and interactions (i.e. 
new physics) are then required. 

In the present work we modify the AF model with an aim to accommodate the above 
desirable features with a parametric relation given by Eq. (6). [We note that deviations 
from tri-bimaximal mixing have already been considered in the AF model, using higher- 
dimensional operators. We take the view of starting with the tri-bimaximal form of the 
neutrino mass matrix given by Eq. (2), however it may arise, and then modifying it with a 
particular radiative mechanism.] The model then predicts new values of the mixing angles 
6i2 and ^13 while bringing the former closer to the central value of the experimental data 
and lifting the latter to nonvanishing values which are within the accessible limits of planned 
experiments [Hj. Because of its structure, any nonvanishing value of ^13 in this model is 
found to be constrained by the solar and atmospheric mixing angles, hence we predict only 
small values of the CHOOZ angle in the end. 

To generate the desirable new radiative contributions to M,^, we introduce three singlet 
charged scalars = 1)2,3) transforming as an 744-triplet in the nonsupersymmetric ver- 
sion of the AF model with two Higgs doublets hu and hd- The Lagrangian of the present 
model has three parts, 

C = Caf + Ci + £2. (7) 

Here Caf is already given in Eq. (1) and in Ref. p], and Ci is the additional contribution 
of the xt scalars that respects A4 symmetry. The term C2 is introduced to break the ^4 
symmetry softly and in conjunction with Ci, it gives rise to new radiative contributions 
known often as the Zee mechanism P], as depicted in Fig. 1. Explicitly Ci and £2 are given 
by 

£1 = / (L L C (3 X 3 X 3) 

^2 = ci2hliT2hd{xt + xt + xt)- (9) 

It is to be noted that L in Ci denotes lepton doublets. The neutrino mass matrix comes out 
as 

fa + 2d/3 -rf/3 -d/3-e \ 

-d/3 2d/3 a-d/3 + e , (10) 

V-rf/3-e a-d/3 + e 2d/3 ) 
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where the a and d terms are obtained in the same way as in the AF model due to higher 
dimensional operators. The e terms are the additional contributions from the one-loop 
radiative diagram as shown in Fig. 1, 

e = fmf-^F{ml,ml), (11) 

Vd 

with the definition 



1 



<h^,> 



V, 



e,fj, 



— « X <- 



Figure 1: One loop radiative i^e^^ - mass due to charged Higgs exchange. 

Apart from the presence of the scalar singlets xt ^ ^^e effective theory below the seesaw 
scale in this case is analogous to the nonsupersymmetric two-Higgs doublet model (2HDM). 
The masses m^^ are degenerate and have been generically represented as m~^. Furthermore, 
the interactions written in Eqs. (8) and (9) can generate corrections to all the off-diagonal 
entries of M^,, but we retain only the dominant terms proportional to m^. Diagonalizing the 
neutrino mass matrix given in Eq. (10) with the assumption that e is small, we obtain the 
three mass eigenvalues as 

mi = a + d + e, m2 = a, = d — a — e, (13) 

where the parameters are, in general, complex. Model predictions for neutrino oscillations 
with Ui3 = have been discussed for complex values of a and d in Ref . jU E] . For the sake 
of simplicity and economy of parameters, we discuss model predictions by treating all the 
three parameters to be real and then, more generally, by treating only d as complex. 

(A) Real parameters 

In this case by treating all the three parameters as real and by solving the eigenvalue equation 
the following expressions are obtained in the leading approximation with |e| << |a|, \d\, 

sin 6^12 = ^ + 5i, sm923 = -{A= + S2), sin6'i3 = 53, (14) 
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where 
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We find that it is possible to fit Amg = m| — and Am^^j^^ — m\—m\ if the two parameters 
are related as 

d — —K a, (16) 
where k is a positive rational number. Then mi — (1 — ft;)a + e, m2 = a, ma = — (l + «;)a — e, 

Aml,^^{K^ + 2K)a\ (17) 



Am| = [(2 - «)/(2 + «)]Am^tm + 2e(« - l)i 
where we have used Eq. (17) to determine the parameter a, 
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In order to estimate the model predictions, we now express e, mass eigenvalues, and mixing 
angles in terms of Am^^j^^, Am^, and the positive number 



((^2 + 2K)Am 
2(fi:- 1) 



2 \l/2 



Amg 
Am^t^ 



K 



2 + K 



mi 



-(«:-ll 



Am^ 



atm 



sin ^12 = — ;= 



k,^ + 2k 
sin 9i3 = 

K + 2 



+ e, m2 



V + 2k 
«;(«; + 3) 



, rn3 = -(«; + l)i 



AmE 



atm 



K2 + 2K 



(20) 
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It is evident from Eqs. (20) to (22) that e and corrections to the mixing angles depend 
upon n apart from the experimentally determined quantities like Amg, Am^^nn ^^^d the 
ratio, R = Am^/Am^^j^. While a positive e would predict sin^^is > 0, sin ^12 < l/\/3, and 
tan^ 6*23 > 1, a negative value would give sin ^13 < 0, sin 6*12 > l/\/3, and tan^ 6*23 < 1. Since 
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the tri-bimaximal prediction corresponding to sin Ou = 1/ ^/S is just on the border hne of 
the maximal value allowed by the recent global analysis [7], the negative values of e and 
sin which shift sin Ou further away are strongly disfavored. We thus search for small and 
positive values of e to predict the masses and mixing angles. 

We note that the number k is not arbitrary. It is clear from Eq. (18) that the smallness 
of AttIq compared to Aml^^ requires k ~ 2. For larger values of k > 2.5 the leading term 
dominance condition, |e| << \a\, \d\, breaks down. For values of 2.2 < n < 2.5, the solar 
mixing angle prediction falls below the present 99% confidence limit and results in 612 < 30°. 
Thus we use the most plausible value n = 2 (i.e. Eq. (6) with = and corresponding to 
the symmetry limit \mi \ = \m2\) and all the parameters in Eqs. (20) to (22) are determined 
in terms of Amg, Aml^^ and their ratio R with 



V2v ^ml^mR, mi = -^-j=^ Amli^ + V2^J Aml^^R, 

3 



^2 = Aml^^, 1713 = -^-j=^Aml^^- V2^Amli^R, (23) 



1 2^6 



sin ^13 = — -^R, sin ^12 = sin ^13, 



4\/2 

tan2^23 = l + -^sin^i3. (24) 

5 

Using the allowed range for Am| = (7.2-8.9) x lO^^ eV^, Am^^^ = (1.7-3.3) x lO^^ eV^, 
we find R = (2.2 — 5.2) x 10"^, and for k = 2, we obtain 

mi = -0.015 eV, ms = 0.017 eV, mg = -0.055 eV. (25) 

Thus the mass eigenvalues are normally ordered jl]. While there is hierarchy between mi^2 
and the masses mi and m2 are nearly quasi-degenerate. The kinematical neutrino mass 
|myj and the effective neutrino mass |mee| contributing to neutrinoless double beta decay 
are also small, 

|m,^J ~ |mi,2| — 0.016 eV, |meo| ^ 0.01 eV, (26) 

which are beyond the detection limits of planned experiments in near future fOlE]- The 
predictions for mixing angles are 

^12 = 31.13° -33.5°, ^13 = 3.5° - 1.5°, ^23 = 45.5° - 46°. (27) 

In Eq. (27), the smaller (larger) value of ^13 is correlated with larger (smaller) value of 6'i2. 
Although still larger values of 6'i3 even closer to the CHOOZ upper limit are permitted by 
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the model they are correlated with smaller values of 612 < 30° and hence are ruled out even 
at 99% confidence level. For example with k = 2.25 we obtain ^13 = 6°, but 612 = 29.2° 
which is below the range allowed at 99% level and hence ruled out. Thus the prediction of 
the angle up to ^^13 ~ 4° is quite natural in this model. The value of z/^ — Ur mixing angle is 
also found to increase slightly beyond the tri-bimaximal prediction. 



(B) One complex and two real parameters 

In this case we treat a and e to be real but d complex with its phase 



arg((i)). In 



order to maintain the experimentally observed smallness of ArriQ compared to Aml^^ we 
use the relation |1] 

|d| = -2acos0. (28) 

Then in the leading approximation, 

|mip = — 2ea(2 cos^ — 1), |^«2p = If^P, 



;i + 8cos2 0)|ap + 2ea(2 cos 



(29) 
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Eq. (28) suggests that lies in the second quadrant for positive values of a. Eqs. (28) to 
(30) give 

'AmU{2V2 cos , 



4cos^ 



(31) 



1 2 cos^ — 1 1 

Thus, the masses \mi\, \m2\, Im^l, and the parameters \d\, \a\ and e are expressed in terms 



eV^ and Am| 



Ami, ^ phase angle 



For example with = 180°, Aml^^ = 2.5 x 10 ^ 



3 X 10~^ eV^, R = 3.2 x 10~^ and we obtain the same values of masses as in 
Eqs. (25) and (26) with normal ordering. With the general expressions for mass eigenvalues 
given in Eq. (11) with complex d, we solve the eigenvalue equation and use Eqs. (28) to (31) 
to derive the following expressions involving the mixing angles. 
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In the suitable limit of | cos(/)| 1, Eqs. (32) to (34) go over , as they should, to expressions 
given in Eq. (22) for the real case with n = 2. We find that interesting solutions bringing 
down the solar neutrino mixing from the tri-bimaximal limit with 612 < 35.3" while increasing 
I sin 1 substantially from its zero limit are possible if the phase of the complex parameter 
d is in the second quadrant. While cos0 = — 1 gives predictions on mixing angles as in 
Eq. (27), Eqs. (32) to (34) can provide substantially different values of mixings for certain 
other values of the parameter, 

cos0 = -0.575: ^12 = 33.5" - 31.2°, ^23 = 45.5° - 46°, ^13 = 1.7° - 3.7°; 

cos (f) = -0.643 : 6^ = 31.5° - 27°, 623 = 45.7° - 46.6°, ^13 = 4° - 8.5°. 

Thus, the prediction for CHOOZ angle could be larger than 4° if 612 < 31.5°. For example 
^13 = 5° would require 612 = 30° which is already near the 2.5cr limit of global analysis. 

The values of light neutrino masses and mixing angles obtained by matching the ex- 
perimentally observed values of AtUq and Am'^^^ are not likely to change significantly by 
radiative corrections through renormalization-group (RG) effects jT2I. This is due to the fact 
that the mass eigenvalues are small |mj| ~ 0(10"^) eV. Further mi and m2 have opposite 
signs and that prevents significant change of sin 612 by RG effects. Although mi and have 
the same sign, they are not so close to produce significant changes in the predicted values 
of sin ^13. In the 2HDM, the radiative corrections are expected to be further reduced in the 
region of small values of tan/3 = Vu/vd ^ 0(1). limit of planned reactor and long baseline 
neutrino experiments |H]. The prediction of solar neutrino mixing angle few degrees below 
the tri-bimaximal value could be tested by precision experiments in near future. Prediction 
of ^13 > 5° is possible if the solar neutrino mixing angle is below the 2.5cr limit of the cur- 
rent global analysis. We find that the present model successfully explains the existing solar, 
atmospheric and CHOOZ experimental results including those from KamLAND and SNO. 
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